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A number of papers over the past eight years have claimed to solve the fractional Schrodinger 
equation for systems ranging from the one-dimensional infinite square well to the Coulomb potential 
to one-dimensional scattering with a rectangular barrier. However, some of the claimed solutions 
ignore the fact that the fractional diffusion operator is inherently nonlocal, preventing the frac- 
tional Schrodinger equation from being solved in the usual piecewise fashion. We focus on the 
one-dimensional infinite square well and show that the purported groundstate, which is based on 
a piecewise approach, is definitely not a solution of the fractional Schrodinger equation for general 
QQ ' fractional parameters a. On a more positive note, we present a solution to the fractional Schrodinger 

' equation for the one- dimensional harmonic oscillator with a — 1. 
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I. INTRODUCTION 



A wide variety of stochastic processes are more general than the familiar Brownian motion, but presumably can 
still be described by modifying the diffusion equation using a fractional Laplacian operator [l|, Q . Such "fractional 
diffusion" is now a large and active field, and a number of books have been written on the mathematics and physics 
of fractional diffusion operators ^i, J,, 5]. In 2000, Laskin introduced the fractional Schrodinger equation, in which 
the normal Schrodinger equation is modified in analogy with fractional diffusion [1, 0, 01 ■ Laskin claimed to exactly 
I '< solve this equation in the case of the one-dimensional infinite square well Q. A more recent (2006) work claimed 
1-^, ' to find solutions again for the infinite one-dimensional square well (agreeing with Laskin's original solution), and for 
d one-dimensional scattering off of a barrier potential \^ . A 2007 work used a different method of analysis to claim 
^ ■ solutions for the linear, delta function, and Coulomb potentials in one dimension [lol |. Laskin also recently built on 
the same claimed solution to derive properties of the quantum kernel The purpose of this work is to point out 
that of the many purported exact solutions presented in the literature, only the one for the delta function potential 

>is correct. 
. The one-dimensional fractional Schrodinger equation 6] is 

^ : zh^^^ = {-h^^Y'^ V(x, t) + V{x, t) ^{x, t), (1) 



where -Da is a constant, A — /dx^ is the Laplacian, and (— Ti^A) ^ is the quantum Riesz fractional derivative: 
Here, 4>{p, t) is the Fourier transform of the wavefunction, 

r+co 

(t>{p,t) = / dx V'(a;,^)e-'P^/^ (3) 



When a — 2, the quantum Riesz fractional derivative becomes equivalent to an ordinary Laplacian, and we recover 
the ordinary Schrodinger equation. 

We focus on the case where the potential is independent of time, so we are interested in solutions of the following 
equation: 

D„ (-a^A)"^^ V(a;) + V(x) iP{x) = EiP{x). (4) 
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The fractional diffusion operator is a nonlocal operator except when a = 0, 2, 4, . . .. This means that {-h''A)'"'\p{x) 
depends not just on ip{y) for y near x, but on tp{y) for all y. This nonlocality, in turn, means that when solving 
Eq. ([4]), the form of the wavefunction in a given region depends not just on the potential in that region, but on the 
potential everywhere. Because of this, for a piecewise-defined potential, we cannot follow the normal strategy of solving 
separately for the wavefunction in each piecewise region, and then using conditions of continuity and differentiability 
to match up the solutions. However, this is precisely the strategy used in the papers cited above, and the solutions 
obtained in those papers are thus invalid. We illustrate the problem by looking in some detail at the case of the 
one-dimensional infinite square well in Section II. The problems with the purported solutions for other potentials 
are similar, and are discussed in Section III. Section IV presents an exact solution for the one-dimensional fractional 
harmonic oscillator with a = 1, followed by conclusions in Section V. 



II. INFINITE ONE-DIMENSIONAL SQUARE WELL 



Consider Eq. ^ in the limit of the potential becoming an infinite square well 

V(r\ - / ° if \x\<a 
~ 1 oo if \x\ > a. 



(5) 



We first note that for the case of free space, where the potential V is zero everywhere, it is easy to see that plane 
waves are eigenf unctions of the quantum fractional Hamiltonian: 



Jpx/h 



(6) 



However, Eq. ([6]) is only valid if the function operated on is e^^^^^ everywhere] it is not a local equation that can 
be applied just in a restricted region. Because the quantum Riesz fractional derivative is a nonlocal operator, the 
wavefunction in the well knows about the wavefunction and potential outside of the well. Previous works looking 
at the one-dimensional infinite square well incorrectly applied Eq. ^ only inside the well, and concluded that the 
solution inside the well would be a simple linear superposition of left- and right-moving plane waves of the same 
energy [l,!!]. 

Although these papers use an invalid assumption, could their end results be correct nevertheless? Both papers 
claim that the solutions for the one-dimensional square well are the same for the fractional case as for the standard 
non-fractional case, only with modified energies. So, they obtain for the ground state 



•00(2;) 



The Fourier transform of this is 



Acos(ff) 



+ 00 



for \x\ < a 
otherwise. 



(7) 



dx e-'P''/''Mx), 

Anh^ cos (ap/h) 

_ (7rn/2a)2' 



From (j^oip) we can calculate the fractional Riesz derivative: 



-h A) ^Jo{x) = — 



dp 



— cos I -7rp 



f'KpX\ 

\2^) 



(8) 



(9) 



We see here how the nonlocality manifests itself in the mathematics. If we only looked at the wavefunction inside 
the square well, then ^"0(2;) would appear to consist of plane waves of just two wavevectors, ±7r/(2a). However, in 
reality, ipoix) is outside the well, making it a wave packet, rather than just a combination of two plane waves, and 
so it contains a continuous range of wavevectors, as seen in Eq. The fractional Riesz derivative thus sees all these 
wavevectors. 

Now we shall show that ipoix) is not a solution of the infinite square well via a proof by contradiction. First, assume 

that ipo{x) is a solution of the fractional Schrodinger equation. Then the fractional Riesz derivative (— ?i^A)" ^ ipo{x) 
must be proportional to ^'0(2;) on the open interval, |a;| < a, where V{x) = 0. Since, ipo{x) is continuous and ipoia) = 0, 
this implies that the limit a; ^ a~ of © should also vanish. However, this condition is not equivalent to ^ vanishing 
at a: = a because the Hamiltonian includes the potential and so we cannot rely on continuity of 4'oix) at x = a. 
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A Fourier transform such as ([9]) is continuous if the integral is absolutely convergent. The integrand in ([9|) is 
bounded by for small p and by (1 + e)p"^^ for large p. Therefore, ^ is indeed a continuous function for all x for 
— 1 < a < 1. Thus, for —1 < a < 1, we can take the limit a; — > ([9]) by setting x — a, and if ^0(2^) is a solution, 
this should give zero: 



/(a) ^ dp cos^ Q^p^ = 0. (10) 

Taking the derivative with respect to a, we see 

The integrand in Eq. (jlip is everywhere positive, so df /da > 0, and wc cannot have /(a) ~ for all a. The ground 
state (O claimed in Refs. @ and thus cannot be an solution of the fractional Schrodinger equation for all a. It 
can only be a solution once in the interval — 1 < a < 1 — namely when a — 0. 

The above argument does not hold for 1 < \a\ and in fact for some values of a, ([9|) is not continuous at a; = a. 
However, a related argument to one presented above shows that this "00 cannot be a solution at least for I < a < 2. 
For a = 2, on the other hand, ipoi^) actually is a solution. Indeed it is a solution whenever the fractional Riesz 
derivative is an ordinary derivative — that is, for a = 0, 2, 4, . . .. 

It may seem counterintuitive that Eq. ([7]) is not the correct ground state. The standard (a = 2) Schrodinger 
equation for an infinite potential well is equivalent to the Schrodinger equation on an interval with the Dirichlet 
boundary conditions ip{—a) "(/-"(o) = 0. By raising that Hamiltonian to the power a/2 we get a plausible fractional 
Laplacian and Eq. ([7|) is indeed a solution. However, this is not the Riesz fractional derivative. In other words, the 
fractional Schrodinger equation for an infinite potential well is not equivalent to the fractional Schrodinger equation 
on an interval. 

At this point, we do not know what the true solutions are for values of a other than 0, 2, 4, . . .. In Ref. flj], Zoia 
et al. find numerical solutions for the ground state. The solutions depend on a and differ from the simple sine wave 
solution in Eq. ([7]). 



III. OTHER SYSTEMS 



While we have only discussed the infinite one-dimensional square well in detail, the comments here equally invalidate 
the other claimed solutions of the fractional Schrodinger equation. For example, in Section II of Ref. [iq|, the linear 
potential, 

. { Fx a x>G , . 

yi^) = ^ ^ if < ' (12) 



00 



is studied. The authors of Ref. ]W{ treat this equation in a piecewise approach by solving the equation for the potential 
V{x) = Fx and applying a boundary condition at a; = 0. This is invalid for the same reasons stated above for the 
square well potential. Similar comments apply to the analysis of the Coulomb potential in section IV of that same 
paper. 

Our comments, however, do not invalidate the analysis of the delta function potential in Section HI of Ref. p^ . 
which did not implement a piecewise ap pro ach, but instead worked with the Fourier transform of the delta function 
potential. However, the authors of Ref. [101 ^^i^ ^o note that the bound state for the delta function potential is valid 
only for a > 1. For a < 1, there is no bound state since the integral in Eq. 33 of Ref. [lO] diverges. In more recent 
work, Dong and Xu attempt to solve the same problem again, using a piecewise approach. They then compare 
this to their initial correct solution and derive an incorrect identity for the H-function. 



IV. THE FRACTIONAL HARMONIC OSCILLATOR 



Consider the fractional Schrodinger equation with the potential 



V{x) = -kx^. 



(13) 
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Fourier transforming Eq. ([4]) gives 

hh'^ = {D^\pr^Emp). (14) 

In momentum space, the equation maps to the ordinary Schrodinger equation with a positive a power law potential 
and k — 1/m. In other words, the kinetic and potential energies have reversed roles. 



A. WKB approximation 

Given the mapping to ordinary quantum mechanics, in the limit h 0, one can use the WKB approximation in 
momentum space to approximate the energy eigenvalues, with p replacing x. The quantization condition in momentum 
space is 

X{p)dp = (n + -)Trh, ne{0}UZ+ , (15) 
pi 2 

where X{p) = ^ |(£' — Da\p\°'), and pi and p2 are the classical turning points, i.e. pi,2 = ±(i?/Da)^/". The above 
condition leads to 

/ (n + l)av^^/fc(J.)V"^(f + l) ^^ 

"""^l 2V2r(i)r(i + i) ) ■ ^''^ 

This agrees with Laskin's more general WKB result [ij for an arbitrary power-law potential. However, by the 
argument we have just given, this special case is better justified than Laskin's general claim. It is unclear whether 
there is any reason to believe the WKB approximation to be vaild for systems other than the harmonic oscillator, since 
for other potentials, the Fourier transform of the fractional Schrodinger equation will not be an ordinary Schrodinger 
equation. 



B. An exact solution for a = 1 

When a — 1, Eq. (fT4|) becomes 

lkh'^^iD^\p\-E)Hp) . (17) 

Restricting to p > or p < 0, this differential equation is equivalent to the Airy equation (by a rescaling transfor- 
mation). For a normalizable wavefunction, we must have (/)(p = oo) = 0(|5 = — oo) = 0. This condition rules out 
Airy functions of the second kind (Bi(z)) as solutions. Because of the symmetry of the potential, the solutions are 
alternately symmetric or antisymmetric. More precisely. 



0(p) = (sgnp)"Ai(AcH-r„) , (18) 



1/3 

where n = (2Di/ (kh^)^ , and the r„'s are the successive roots of Ai' (for n even), or of Ai (for n odd). The energy 
eigenvalues are 

rk \^/^ 

£;„ = - [^^n'Dlj r„. (19) 

Note that r„ < 0, so > 0. 

Using the asymptotic expansion of the Airy function, we find that the roots are well approximated by 

o / , X X 2/3 

Stt f 1 ^ ^ 



^ ^n+-jj , ne{0}UZ+ (20) 

in the limit of large n. This approximation reproduces the result of the WKB approximation above. This approximate 
formula is off by 8.7% for n = 0, 0.77% for n = 1, and 0.41% for n — 2, and rapidly becomes more accurate for larger 
n. Inserting Eq. (I20p into Eq. p9p thus gives a very good approximate formula for the energies of the simple harmonic 
oscillator for a — 1. 
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V. CONCLUSIONS 



It would be useful to know the correct groundstate of the one-dimensional infinite square well or harmonic oscillator 
for general a, but this is a difficult problem. In Ref. [lB|, Bahuelos et al. needed a lengthy proof merely to show that 
the groundstate solution for the infinite square well in the region (—1, 1) is concave on the interval (—5,-1-5). 

Similar technical issues regarding nonlocality have arisen in the statistical mechanics community as well. For 
example, in Ref. [IB], the authors claimed to analytically determine the mean first passage time for a Levy flight 
with absorbing boundary conditions on the interval [0, 1]. They did so by imposing the standard absorbing boundary 
conditions for the probability density at a; = and x ~ L. However, a subsequent publication pointed out that 
due to the nonlocal nature of the Levy flight, the correct boundary condition is for the probability density to vanish 
for all a; < 0, and for all x> L, rendering the analysis in Ref. 16 ] invalid. 

Finally, one must also ask about possible physical realizations of the fractional Schrodinger equation. In Ref. 
[13], Zoia, Rosso, and Kardar constructed a lattice model whose continuum limit is described by fractional diffusion 
using symmetric Toeplitz matrices. A quantum representation of this model via a mesoscopic network of long range 
connections whose hopping amplitudes are described by the entries of these matrices may be a realization of the 
fractional Schrodinger equation. Further modification could lead to experimental tests of the problems we have 
discussed here. 

JMS would like to thank the hospitality of the Aspen Center for Physics where part of this work was completed 
and support from NFS-DMR-0645373. 
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